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Abstract 

I briefly review the three nonperturbative methods for the treat- 
ment of disordered systems — supersymmetry, replicas and dynamics 
- with a parallel presentation that highlights their connections and 
differences. 



Disordered systems need to be treated with a method that allows to 
perform averages over the sample realisation. There is no universal way to 
do this that can be applied efficiently to all problems. 

For Gaussian systems, the method of supersymmetry is as good as one 
can expect: it involves a minimum of variables, it is elegant and rigorous. 
Although one can still apply it for some non-Gaussian problems, in many 
of the interesting cases - as for example spin-glasses - it only gives limited 
information. 

The replica trick was introduced to tackle such 'complex' problems. It has 
been extensively used and has yielded some of the most innovative solutions 
in disordered systems. It has however the problem that it is very far from 
being controlled, let alone rigorous. This is because the space itself - a vector 
space with noninteger dimension - does not have a general definition other 
than the ansatz itself - or perturbations around it. 

The dynamic method consists of solving exactly the evolution of the sys- 
tem in contact with a heat bath. If the system reaches equilibrium one 
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recovers all the thermodynamic information. Surprisingly enough, one can 
treat this way all the problems one can solve with replicas. A problem arises, 
however, when equilibrium cannot be achieved: then the long-time out of 
equilibrium regime may be of interest in itself (as in the case of glasses), or it 
may be viewed as an obstacle for exploring the deepest levels in phase-space 
(as for example in optimisation problems). Although the dynamic method 
was initially proposed as a way to obtain equilibrium results, this tendency 
has reverted in the last few years, at least in the field of glasses, where replicas 
are now used mostly to mimic the out of equilibrium dynamics. 

The aim is of this paper is not to make a complete presentation of either 
of the three methods — there are very complete reviews of this || |3|] 
(including some very recent ones M) but rather to put the three methods 
'side by side' so that the connections can be better appreciated. To the best 
of my knowledge this has not been done for supersymmetry, replicas and 
dynamics simultaneously, as the practitioners of each method tend to belong 
to different communities. 



The Problem 



Consider an energy 

Ej = % - y^(A^ 3 - - Jij)siSj ; Ej{h) = Ej - ^ h i s i 



'j 



where Si (i = 1,...,N) are real variables, and is a random matrix. We 
take A with negative imaginary part. This energy can be used to calculate 
the averaged Green function: 



G(X) = Tr[AI - J] -1 



(2) 



from which one obtains the eigenvalue distribution. (Here and in what follows 
the overline denotes averages over the disorder J). This is done by defining 
the partition function 



Zj(h) = J ds e-^W 



(3) 



and computing: 



G(A) 



d 2 



-* T J2w lnZj{h) 



hi=0 



ds s 2 k e-^M 



Zj(h) 



ds s 2 k e-^W 



Zj(h) 



hi=0 



(4) 



hi=0 
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(We include the constant normalisations in the differential: d = dj y/j2irN)). 

The third expression in @ is the correct (quenched) average, in general 
different from the last one, the annealed average. The problem is that in 
order to compute the average over the Jij, we need to express 1/Zj in (^) in 
a tractable (i.e. exponential) form. Three methods to do so are: 

• Supersymmetry: we can take advantage of the Gaussian nature of the 
partition function to write, in terms of two sets of Grassmann rji and 
77* and a set of ordinary variables af 
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z] 



dr] drfda e -f i £i^Wy)W'w+^) ( 5 ) 



so that we get: 



G(X) = /Zj dsd(T j dr > dri * s k e~^ iEi i (XSir 



hi=0 

(6) 



Replicas: we replicate n times each variable Sj — > sf and compute: 

Zy 1 = J H n a - = \ds a e-f^lE^-ASy).?-? (7) 

The calculation proceeds for every integer n, and finally we somehow 
take the limit: 

Z -j 1 = lim Zy 1 (8) 

which should in principle be shown to be the correct analytic continu- 
ation over n. We hence have: 



G{\) ~ n _ [ Ka=ids a (4 X) ) 2 e-tS-iE^-AawK--? (9) 
k J 

where we have chosen to take the expectation value of the first replica, 
although clearly any other replica will do. 

Dynamics: The dynamic method || || consists of calculating the aver- 
age in (f|) by considering the solution of the Langevin equation: 

7s i = -P-Q±+/3h i + p i (10) 
where the pi are independent Gaussian white noises with variance = 27. 
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The energy might be complex: this poses no problem (at least for linear 
systems f7j). Starting from t = 0, we are guaranteed that at long times 

(A(s)) = Vm(A{a{t ))) p (11) 

where (•) denotes thermodynamic average and («) p average over the 
process, i.e. over the noise realisation. We obtain an expression for the 
average Green function (f|) as: 



G(X) = i lim V 



d(s k (t )) f 



to^oo^-^ dhi 



(12) 

h=0 



In practice, one calculates the dynamics averaged over both thermal 
noise and disorder and in the large N limit, as we shall see below. 

The problem of treating the denominator is not exclusive of Gaussian 
systems, it appears whenever we wish to obtain the correct quenched averages 
over disorder. For example, the energy ([]]) can be modified to obtain the 
standard spin-glass model: 

E n J l = i E J (h)+mJ2s 2 l +gJ2 s t (13) 

i i 

and we may wish to calculate averages of any observable A(s). As soon as 
g > the system becomes as complicated as can be, with all the subtelties 
of spin-glasses. Once we abandon the Gaussian world, the three methods 
encounter difficulties: 

• Supersymmetry: there is no obvious way to write 1/Zj in general as an 
integral over an exponential. This does not mean that the supersym- 
metry method is entirely inapplicable for non-Gaussian systems: even 
though when the energy is not quadratic this method does not give the 
Boltzmann-Gibbs measure, it can still be useful in some cases, as we 
shall see below. 

• Replicas: In contrast to (H), expressions (f7|) and (g) are formally valid 
for non-quadratic energies. Thus, the replica trick has been applied 
successfully to the study of many complex systems, spin-glasses being 
the main example. The expectation values of an observable A can in 
general be written as: 



(A) ~ n _ / Tl n a=1 ds a A(s^) e -^.^(«°) (14) 
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From the point of view of making the results rigorous (or even reliable), 
there is the following difficulty: a closed analytic expression in terms 
of n can be obtained in some limit, typically large N. This poses the 
problem that the limits N — > oo and n —>■ may not commute - and 
indeed in most interesting cases they do not. In those cases we have to 
consider the assumed infinite- N continuation valid around n = as a 
guess (see however Ref. fit 



The dynamic expression (O) shares with the replica treatment the 
advantage of being equally valid for linear or nonlinear problems. There 
is however a problem also here: ([11]) holds to the extent that we make 
t a — > oo before any other limit, in particular iV — > oo. Again, in 
many interesting (nonlinear) problems these limits do not commute: 
in physical terms this means that an infinite system is not able to 
equilibrate at finite times M. This is indeed the physical situation 
one wishes to reproduce in glassy systems. However, one may still 
be interested in knowing what happens in times that diverge with the 
system size, and in particular to reproduce the equilibrium situation - 
even if it might be unreachable in a realistic situation Q. To do this, 
the N —> oo solutions must be supplemented with activated, 'instanton' 



solutions [11]: this problem has not yet been solved in general. 



Dynamics is a generalisation of supersymmetry. 



Let us see that the supersymmetry method is a 'time-less' version of 
dynamics (PHI). We compute the solutions of the stochastic equation: 



dE 

= -P-^ + Ph + Pi 

OSi 



(15) 



There is no time-dependence, and the pi are Gaussian variables of variance 
27. If Ej{h) is quadratic the system flloD has a single solution 



(16) 



Denoting (A(s)) the average of A evaluated over the (p-dependent) solutions, 
we have : 



d(s k ) f 
dh k 



(17) 



h=0 
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to be compared with (12]). To see that this gives back the supersymmetry 



method, let us write, for the Gaussian case 



J ds s k UiS J2( X6 v ~ J v) s i + P h i + P^j det[i/3(AI - J)]j 

(18) 

where the determinant guarantees that the solution for every realisation of p 
is counted with the same weight. Exponentiating the delta function as usual 



12 : 



(sk)p — (^j dsdsdrjdrj* s k x 

exp <^ ^ -i/3[\5ij - Jij}(isiSj + 77*%) + i ^ Si(f3hi + p { ) > 

I ij i ) 

= J dssd-qdrj* s k x 

exp < ^2 - Jij](isi-Sj + r]*Vj) + i ^ ~~ 7 ^2 ^ \ 

(. ij i i ) 

(19) 

which, using fllTD yields: 

G{\) = i/3 22 / dsdsdrjdrj* s k s k x 
k J 

exp < Y -iPfiSij ~ JijWsiSj + rj*rjj) - 7^ s 

I ij i 



(20) 



This is an implementation of supersymmetry like @, with two ordinary (s, s) 
and two Grassmann (77, r?*) sets of variables. For 7 — > it can be taken to 
the form (|6|) by a rotation in the (s, s). 

The conclusion we draw from this exercise is that: i) Supersymmetry is 
just 'dynamics without time', which strongly suggests that any problem solv- 
able with the former is solvable with the latter method, ii) Supersymmetry 
can be extended to treat certain nonlinear problems, as we shall now show. 

Supersymmetry for nonlinear problems. 

Equation (|1^) is not restricted to linear energy functions. If (|5|) is non- 
linear, but still has one solution, it can be used to calculate the expectation 
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value of any function A(s) in its root. The generalisation of eqs. (|T8D and 



(0) is: 

(A(s)) p 



,1s A U,S { -P^- +Pi ) (] ^ 



d 2 Ej 

ds k dsi 



dsdsdrjdt)* A x 



dE 



exp 



ds,;ds 



-v., 



(21) 



This way of imposing a solution has its origin in the path-integral treatment 



of gauge theories | 13| , where the fermions are called 'ghosts' [fL4 |. 

In many cases of interest the equation ( p~5| ) has many solutions for some 
realisations of p. If we wish to add the values of the observable in every 
solution we should take the absolute value of the determinant in (pip. In 
particular, we need to do this if we wish to calculate the average number 
of solutions. Writing this absolute value as an exponential is possible P~5|] , 
although it involves introducing new fields. 

An interesting situation we shall consider here and in what follows is when 
we do not take the absolute value. Each solution is then added with the sign 
of the determinant of the matrix of second derivatives [Q. In particular: 



(i)= E 



sign 



(22) 



solutions 



which is an invariant only dependent on the topology of the space of the s, 
and independent of the energy function Ej 
s forming a flat space and Ej(s) — > oo as 



For the usual case of the 
\s\ — > oo the invariant is one. 
In cases in which there are many solutions, the method does not select the 
lowest ones, but averages flatly (apart from the sign of the Hessian) over all 



solutions [fl7 |: it is in this sense that supersymmetry fails. 

In any case, as mentioned above, one is not calculating the Gibbs measure, 
but just values over local minima and saddles. There are however interesting 
nonlinear problems having a finite number of solutions for which there is no 
reason to abandon supersymmetry. 

One of the most interesting applications involving non-gaussian prob- 
lems are the quantum systems. A note on terminology is necessary for what 
follows: In quantum systems, we can distinguish two ways in which non- 
linearity may appear: in the wavefunction and/or in the Hamiltonian. In 
the former case, one has a nonlinear Schroedinger equation, containing for 
example terms cubic in the wavefunction (see Eq. (121) below). In the latter 



7 



case, one generally considers a usual, linear Schroedinger problem, but the 
Hamiltonian contains terms of degree higher than two in the creation and 
destruction operators. It is then the path integral that is non-Gaussian, since 
the action is no longer quadratic. We shall discuss below both cases. 



Functional expression for dynamics. 



We can see more clearly the relation between supersymmetry and dy- 
namics by constructing a functional expression for the equation (TO). We 
use exactly the same procedure as in (]21f) , with now delta-functions and 
Jacobians promoted to functionals of the trajectories. 

(A(s(t ))) p = J DsDsDrjDri* A(s(t )) x 

exp j -tf £ J 'dt St^L + /? £ J dt »tf|J^ 

+7^1 dt^-ihk -si) J (23) 
This functional equation can be viewed either as the de Dominicis-Janssen, 



Martin-Siggia-Rose [|T8|, |T^] functional expression for the Langevin dynam- 
ics - with the determinant exponentiated through ghosts - or as the path- 
integral expression for supersymmetric quantum mechanics [EDI . 

Here we see clearly that by expressing expectation values dynamically the 
problem now becomes, just like in the case of supersymmetry and replicas, the 
computation of an integral of an exponential, albeit a functional one. This is 
the usual starting point for the developments in dynamics - at least within 
the physics literature. 

This is a good place to see how one can calculate with the same method 



the localisation of wavefunctions in a nonlinear Schroedinger problem ]21 

ifin = ~((Pn-l + <Pn+l) + (fin + MVnl^Vn (24) 



where e n is the site disorder. The system of equations fl24j) together with 
its conjugate can have a single solution (this will surely be the case if we fix 
the initial conditions). In fact, ( 21 ) can be viewed as a (noisless) dynamical 



equation for the variables tp n . We can obtain a functional expression as in (|23|) 
for this case introducing complex Lagrange multiplier time-dependent fields 
4> n {t), and Grassmann fields r) n (t), T]*(t). The normalisation is guaranteed, 
even if the action is no longer quadratic in the 4> n {t). 
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Normalisation and symmetries. 



We have three expressions for the expectation of an observable: using su- 
persymmetry ([H]), replicas ( |HD and dynamics (|23|). All three lend themselves 
to averaging over the disorder, and have no uncomfortable normalisations. 
Indeed, the three expressions yield 

(1) = 1 (25) 
but for apparently different reasons: 



Within the supersymmetric formalism ( p5|) arises because around each 
solution the Grassmann and the ordinary variables conspire, just as 
in the Gaussian case, to give ±1 (the sign of the determinant of the 
Hessian). Even when there are many solutions, these signs add up to 
one because of topological constraints 



Now, even if we did not know where the function fl21|) came from, we 
could still see that the expectation value (1) does not depend on Ej 
using the fact that the exponent has the two supersymmetries (which 
indeed give the name to the approach): 

5si = rji ; 5r]* { = isi 

5si = rj* ; Srji = isi (26) 



Within the replica formalism (^) just expresses the fact that we have 
an integral to the n power, and we let n — > 0. Again, if we did not 
know where (1T3) came from, we could show that (1) = 1 using the fact 



that the exponent is symmetric with respect to replica permutations. 

In the causal dynamic treatment starting from an initial condition and 
letting the endpoint free, (|25"D is just a statement of probability con- 
servation f35fl . Also in this case we can see directly from the action 
that (1) = 1, for reasons of symmetry One has the following two 



supersymmetries, which are the generalisation of (|26|) to the case 'with 
time ': 

5si = r]i ($77% = isi 

5si = rj* 5r)i = isi - si 5si = -irj* (27) 

which, together with time-translation invariance, constitute the full group of 
symmetry. 
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A unifying notation. 

We have seen that the methods of supersymmetry and dynamics (itself 
also possessing a supersymmetry) are closely connected. In fact, we can 
uncover more algebraic correspondences between the three approaches by 
using a suitable notation [|T^, |23J. This can be done by introducing two 
anticommuting Grassmann variables , 0: 

[0 ; 0j + = e 2 = 2 = (28) 

The integrals over these variables are defined as: 

J 1(20 = J ld9 = J 9d9 = J 9d9 = 1 (29) 

We can encode the Sj, rji, rj* and s, in a single superfield: 

$. =s . + g ri . + v *g + §.g ( 30 ) 

Using Eqs. (P8|)-(|2"9D and ( |3"U] ) one obtains, in terms of the superfields $j 

(A) = JY[D[^]A exp J da ^ ^(a)A(a, a')$i(a') - PEj($ (a)) 

(31) 

where we have denoted a = (0, 0), cfo = <i0 <i0 and 

A = A SUSY (a,a') = 2 1 (32) 

independent of a, a'. 

The dynamics can be encoded in an expression formally identical to (|31"D, 
but now the field dependencies and integration variables include time: a = 
(9, 6, t),a'= (9', 0', *'), da = d9 d9 dt and 

A = A. Dyn = 2j6{t - t') + j6'{t - t') (9 - 9') (0 + 0') (33) 

Finally, the replica expression is again formally ([H]), but with the identifica- 
tion: 

/n 
da ( 34 ) 

(The correspondence between supersymmetry and replicas can be made to 
hold even for 7 7^ by using a term A p = 7 which does not affect the 
final result.) 
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In particular, the expectation values (A) associated with the calculation 
of the Green function (j9|),(|12) and (20) can be written in terms of A which, 



in a notation that highlights the analogies, reads: 

A = i^Y^ / dada ' <5>a{a)0{a')<5> a {a') (35) 

with the identifications: 

0(a) = 8 aA ; 0(a) = 5(6)8(6) ; 0(a) = S(t - 1 )8(6)6(6) (36) 

for the replica, the supersymmetry and the dynamic cases, respectively. We 
see that the expressions are analogous to one another. 

The important point about expression ([31]) is that, apart from the first 
term in the exponent, it has the same form as the partition function. This 
unified notation is useful as a book-keeping device when we have a diagram- 



matic expansion p3[ [24j| , because diagrams on the three approaches have the 
same form. Internal lines involve integrations over the superspace/replica 
variable, and the effect of each method is the same due to relations like: 



1 da = (37) 

valid in all three approaches (in the replica approach as n — > 0). 
The correspondence at work. 



One-point functions of random matrices. 

We now work out the example of the one-point function for the Gaus- 
sian orthogonal ensemble in parallel with replicas and supersymmetry. (The 
problem has also been attacked with dynamics f25| , but we will not review 



this here). The object is not to discuss how both methods can be used in 
this case (this has been done in detail long time ago |41|, |26|), but rather to 
show how the equality of results follows from the formal correspondence. 

We use the functional expression fl3~I|) with the energy given by (|I|) , where 
the Jij are random Gaussian variables of variance iV -1 / 2 . Averaging over the 
J, and expressing everything in terms of the order parameter: 



Q(^ a ') = -^X>i(a)$i(a') (3* 
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we get, after a few standard steps (which can be borrowed either from the 
supersymmetry or from the replica literature): 



G(X) = J dada'0(a')(Q(a,a')) 
(Q{a,a')) = [ D[Q] Q(a,a') x 



exp j-y Tr \n[AS + i/3X + f3 2 Q] + j [3 2 Tr Q 2 

' (39) 

Here we have used (p5|). The square and log functions are to be understood as 
functions of Q considered as an operator (i.e. Q 2 (a,7) = J da'Q(a,a')Q(a', r y), 
etc) and Tr Q = J daQ(a,a). The delta function is either the Kroenecker 
function (in replica space) or the superspace delta S — (9 — 0'){6 — 6'). 
Expression (|39D can be evaluated by saddle point integration. 

Q- 1 = (A + i{3\)5 + (3 2 Q (40) 

We can now propose for the saddle point value the most general (replica and 
super) symmetric form for Q: 

Q(a,a') = qd + q (41) 

First note that under operator powers and traces (fll]) behaves exactly in 
the same way whether we interpret it as being a replica matrix (n —>■ 0) or 
as a function of two superspace variables. The saddle point equation then 
becomes: 

1 = if3\q + l3 2 q 2 

= 1 q + i(3\ q + 2(3 2 q (42) 

Using (|35|) we have: 

G{\) =i(3 J dada' Q(a, a')0(a') = i(3q (43) 

and this yields the semicircle law in the usual way. 

The point worth noting here is that there is a close algebraic relation 
between the replica and the supersymmetric approaches. Indeed, as we shall 
stress below, all three approaches are essentially isomorphic when restricted 
to a symmetric ansatz. 
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Quantum systems with interactions. 



As a second example, let us briefly see how dynamics can be used as 
an alternative to replicas in an interacting quantum system. Consider the 
system of interacting bosons in a random potential |27 with imaginary-time 
action: 

S = j d 2 xdr^* (d T - ^-V 2 - n + V{x)^Ji) 

+ I d 2 xd d x'drtfj*(x)ilj{x)u{x-x')tfj*(x')4j{x') (44) 



where u(x — x') is the boson interaction and V(x) is the random potential. 
In order to do the correct averaging over disorder, one can use the replica 
trick, thereby obtaining the averaged action: 



S = J d 2 xdrtjj* a (x,r) (d T - ^- V 2 - fj\ ^ a (x , r) 



d 2 x dr dr' ~v ^* (x, T)i[> a (x, r)if)p(x, r)ipp(x, r') 

+ d 2 xd 2 x'dTipl(x)ip a (x)u(x-x / )ipl(x')ip a (x') (45) 



a — 1, 2, . . . , n is a replica index. 

We can just as well apply a dynamic treatment here. Going back to 
fl44l) , we can consider x and r as the site indices, ip(x,T) and ip*(x,T) as 
the dynamic variables, and consider their Langevin evolution in an extra 
(unphysical) time t: 

dt 5^(x,r) +P[,T,t) 1 j 

This 'stochastic quantisation' strategy can be implemented for fermions as 



well [p8 |. We can obtain an expression that is formally identical to (^) (up 
to a term A as in (j32|) ) , but now interpreting the fields ip as superfields, 
functions of both x, t and the superspace variable a = 9,9, t. Diagrams for 
superfields have the same form as the replica ones, and one can also study 
nonperturbative approximations. 

Let us conclude this section by remarking that for this last case there is 
another (more physical) approach: the treatment of quantum dynamics with 
a thermal bath a la Schwinger-Keldysh (see the first of Refs. [|J). This has 
the advantage of not having to introduce an extra time. 
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Order parameters, symmetry breaking. 



Order parameters can be of vector nature ^(cc), of matrix nature Q(a, a') 
and of higher tensorial character. They may, of course, depend on space. A 
special case arises when one wishes to calculate the two-point correlation 
function of random matrices. One needs to introduce two sets of superfields, 
or of replicas $^(a), &\ 2 \a), and ends up with an order parameter: 



n _(Q {ll) Q (l2) \ 
Q ~ V Q {21) Q {22) ) 



(47) 



where NQ^ (a, a') = } \oc')) for a, b = 1,2. 

The different solutions can be classified according to the manner in which 
the symmetry is broken. 

• Symmetric order parameters appear in the solution of Gaussian one- 
point problems. This corresponds, as we have seen in the previous sec- 
tion, to replica-symmetric/supersymmetric solutions. In the dynamic 
treatment, the fact that correlation functions satisfy supersymmetry 
(P7|) is equivalent to stating that the system is in equilibrium, and sat- 
ifies stationarity as well as the fluctuation- dissipation theorem. The 
dynamics of glassy systems in the high temperature phase is of this 
kind, and can be solved easily in all the cases in which the replica 
trick calculation can also be implemented. (For an explicit presentation 
of the algebraic connection between the two methods, see JOJ). 

• Vector breakings 

Within the replica trick such form of symmetry breaking appears when 
the order parameter is a vector in replica space, and all components 
are not equal For matrix order parameters, vector breakings are 



those such that the vector \1/ defined as: 

»(„> S /**,(„, o>) (48) 

is itself non-symmetric, i.e. dependent on a. The same definition can be 
applied to supersymmetric and dynamic solutions, with the substitution 
of 'replica-symmetry' by 'super-symmetry'. There are several examples 
of such symmetry-breaking fields in the literature: i) vectors in replica 
space [^] were considered in the study of instantons in the random field 
Ising model, their supersymmetric and dynamic counterparts [^TJ have 
closely related properties, ii) Replica matrices with vector type were 
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considered 32 in the computation of saddles in free-energy landscapes, 



and also in for the two-point functions for random matrices. A 
related scheme with matrices is the 'two block model' |3JI], (the first 
attempt at replica symmetry breaking) used to count solutions of a 
spin-glass equations. For this last example there is a supersymmetry- 



breaking ansatz shown to have the same properties [34, 15|, and more 



recently a causality-breaking dynamics |p6 



Matrix breakings: This appear only for two (or more) indexed corre- 
lations. They can be characterised by the fact that although Q(a, a') 
breaks the symmetry, the integral \& (Eq. (|4*8|)) is itself symmetric (in- 
dependent of a). The best known example of matrix breaking is the 
Parisi ansatz in replica space. In the context of dynamics the solu- 
tion of the long-time out of equilibrium evolution of the same systems 
P7| , |3| is of this kind. Both the Parisi ansatz and the dynamic solution 
have been generalised to order parameters of higher tensorial character 



SHI 



Whenever the replica trick is feasable, the dynamic treatment is also 
possible. They do not yield the same answers if the system is not er- 
godic, as one corresponds to the equilibrium situation and the other 
to the nonequilibrium dynamics. Only with the inclusion of all ac- 
tivated (instanton) processes will the dynamic solution reproduce all 
time regimes, and this is not yet available in general O, |36| . 



In several of the cases above, the equality between the solutions within the 
different methods stems from an algebraic correspondence, a generalisation 
of the kind of that we described in the previous section. 



Conclusions 



Having a dictionary that allows to translate developments from one method 
to the other, whenever this is possible, can be useful for several reasons. For 
example, in the field of structural glasses and supercooled liquids, arguably 
the most important theoretical challenge is the inclusion of solutions repre- 
senting the activated processes responsible for the smearing of the purely dy- 
namic transition. Once these solutions are found, one can envisage construct- 
ing formally analogous solutions in replica space, which one might conjecture 
would be responsible for the disapearence of the thermodynamic (Kauzman) 
glass transition, or for a change in its nature. 

From the point of view of mathematical physics, the dynamic method 
seems a promising strategy, since everything that is involved is standard 
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probability theory and analysis |40], [43]]. Indeed, there seems to be no obstacle 
of principle for the rigorous derivation of the solution of out of equilibrium 
glass dynamics []3"7| , f|2, Rl, at least at the mean- field level. 
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